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of states, and includes a huge enhancement for photon frequencies near the charge-density-wave gap 
energy. 

PACS numbers: 71.10.Fd, 71.45.Lr, 78.30.-j 



I. INTRODUCTION 

Inelastic light scattering is a powerful probe of the 
charge fluctuations in a strongly correlated material^. By 
using polarizers on the incident and the reflected light, 
one can examine different symmetry channels for the 
charge excitations, and how easily they can scatter light. 
Using inelastic light scattering, one can learn about the 
symmetry of underlying order, such as the d-wave super- 
conductivity in the high temperature superconductors. 
Here, we will focus on the effects of static charge-density- 
wave (CDW) order on the inelastic light scattering of a 
strongly correlated material. 

The field of inelastic light scattering has been increas- 
ing in interest. When x-rays are used for the light 
source, one can examine resonant inelastic x-ray scatter- 
ing, where both energy and momentum are exchanged 
between the light and the charge excitations of the solid. 
Here, we focus on the zero momentum limit, where only 
energy is exchanged, because we will be using optical 
light. Hence we will be examining resonant effects in elec- 
tronic Raman scattering. The dynamical mean-field the- 
ory (DMFT) approach to this problem was completed a 
few years ago£~— in the normal state. One of the interest- 
ing results from that work was that one could see a joint 
resonance of low-energy features, with higher-energy fea- 
tures when the photon energy was on the order of the 
interaction strength U between the electrons. When one 
has charge-density-wave order, there are two additional 
complications that arise: (i) the density of states (DOS) 
has significant temperature dependence below T c , where 
excitations with energies smaller than the gap energy will 
be depleted as T —> 0, and (ii) the DOS develops sharp, 
singular peaks as T — > that arise at the gap edge. One 
would hence expect the Raman response to have much 
more temperature dependence than what was seen in the 
normal state and to have more striking resonant effects 
because of the sharp peaks which develop due to a pile- 
up of the DOS at the gap edge. Indeed, the nonresonant 
Raman response, in the CDW phase, shows dramatic ef- 



fects due to the singularity in the DOS in some of the 
symmetry channels 16 . 

CDW order is also interesting because there are a num- 
ber of strongly correlated materials that display this be- 
havior. The most prevalent class of such materials are the 
transition metal di- and trichalchogenides, which display 
either quasi one dimensional (NbSes) or quasi two di- 
mensional (TaSe 2 or TbTe 3 ) CDW order££. In addition, 
there are known three-dimensional systems like BaBi03 
and Bai-^K^BiOa which display charge-density-wave or- 
der via nesting on a bipartite lattice at half filling 9 . This 
latter example is particularly relevant to our work, since 
the DMFT is more accurate as the dimensionality in- 
creases. One of the longstanding questions in the field is 
the question of whether the order is driven electronically, 
with a lattice instability following the electronic instabil- 
ity or vice versa. We won't have any direct answers to 
that question in this work, since we are not examining 
time-resolved phenomena, but we will note that exper- 
imental light scattering work has already examined the 
phonon softening phenomena that is associated with the 
lattice distortion^ .. Here we focus on electronic effects, 
which would be the obvious next generation of experi- 
mental probes on these systems. 

We will be varying the photon energy over a wide range 
of different values. We will see the most remarkable res- 
onant effects when the photon energy is equal to the gap 
energy as one might naively expect. For many CDW 
systems, this gap energy is at most a few hundred meV, 
which is much below the optical photon energies. Hence, 
the experimentally most relevant results will rely on ex- 
amining joint resonant effects, like what was observed 
in the normal state in previous calculations. But we 
also will focus some attention on the most dramatic res- 
onant effects under the hope that such CDW systems, 
made from strongly correlated electronic systems, might 
be found in the future, and that they can be studied with 
electronic Raman scattering. 

We use the Falicov-Kimball model in our analysis be- 
cause it is one of the simplest models^ which possesses 
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static CDW ordering and has an exact solution within 
DMFT 12 (for a review see Ref. [HI). In particular, the 
irreducible charge vertex is known exactly, and that is 
needed to examine the charge screening effects. Our 
work also extends recent results on transport, optical 
conductivity, and nonresonant x-ray scattering in CDW 
systems^— to the realm of resonant inelastic light scat- 
tering. A brief report on resonant Raman scattering has 
also appeared as a conference proceedingil. 

The organization of this paper is as follows: in section 
II, we introduce the model and briefly review the dynam- 
ical mean-field theory approach in the ordered phase; in 
section III, we describe the general formalism for inelastic 
light scattering; in section IV, we focus on the detailed 
formulas for the mixed and resonant contributions to Ra- 
man scattering; in section V, we present our numerical 
results and we analyze the Raman scattering response 
for two different cases; and in section VI, we present our 
conclusions. 



II. ORDERED PHASE DYNAMICAL 
MEAN-FIELD THEORY 

Historically, the Falicov-Kimball modelii was intro- 
duced in 1969 to describe metal- insulator transitions in 
rare-earth compounds and transition-metal oxides in- 
volving a simplified two-band model with localized heavy 
electrons and itinerant light electrons which hop between 
sites. The mobile electrons hop to neighboring sites with 
a hopping integral —t and they interact with the local- 
ized particles at the same site with the Coulomb energy 
U. The mobile electron creation (annihilation) opera- 
tor at site i is denoted by d\ (dj and the local electron 
creation (annihilation) operator at site i is fj (/J. We 
perform our calculations at half-filling because, in this 
case, there is an insulating CDW phase at low temper- 
ature for all values of U . The explicit formula for the 
Hamiltonian appears in Eqs. (Q]-[2]). 

An algorithm to determine the (period-two) ordered- 
phase Green functions (within DMFT) was developed by 
Brandt and Mielschi^ shortly after Metzner and Voll- 
hardt introduced the idea of the many-body problem 
simplification in large dimensions- 9 . The CDW order 
parameter displays anomalous behavior at weak cou- 
plin g 20 ' 21 , and higher-period ordered phases are possi- 
ble, and have been examined on the Bethe lattice^. In 
previous worksi^— , the transport properties and nonres- 
onant inelastic light and x-ray scattering were examined 
in the commensurate CDW phase. A detailed descrip- 
tion of the DMFT solution for the CDW phase of the 
Falicov-Kimball model has also appeared in our previous 
paper s 15 ' 16 , so we restrict ourselves to a brief summary 
in order to establish our notation. 

We work on an infinite-dimensional hypercubic lattice 
with nearest neighbor hopping. This lattice is bipartite, 
implying that it can be divided into two sublattices, de- 
noted A and B, with the hopping being nonzero only be- 



tween the different sublattices. In this case, the Falicov- 
Kimball model has particle-hole symmetry, and the non- 
interacting Fermi surface is nested at half filling with 
an ordering wavevector at the zone boundary along the 
diagonal, which implies the CDW order will lie on the 
sublattice structure, with the density of the light and of 
the heavy electrons being uniform on each sublattice, but 
different on the different sublattices. This difference in 
electron filling serves as the order parameter for the CDW 
phase. Keeping this in the mind, we introduce sublattice 
indices into the Falicov-Kimball model Hamiltonian 

H = Y,H? -Y,t<$$J jb , (1) 

ia ijab 

where i and a = A or B are the site and sublattice in- 
dices, respectively, and tfj is the hopping matrix, which is 
nonzero only between different sublattices (tfj A — t BB — 
0). The local part of the Hamiltonian is equal to 

H? = Un? d n? f -rtnt d -f, a f nt f -, (2) 

with the number operators of the itinerant and localized 
electrons given by hid = $ i d i and hif = respec- 
tively. For computational convenience, we have intro- 
duced different chemical potentials for different sublat- 
tices, which allows us to work with a fixed order pa- 
rameter, rather than iterating the DMFT equations to 
determine the order parameter (which is subject to crit- 
ical slowing down near T c ). The system achieves its 
equilibrium state when the chemical potential is uniform 
throughout the lattice (/j, a = \i B and /j, a — A 1 / )• 

The first step of the DMFT approach is to scale the 
hopping matrix element as —t = —t*/2yD (we use t* = 1 
as the unit of energy) and then take the limit of infinite 
dimensions D — > oo.~ The self-energy is then local: 

E?/(w) = £»<%<U, (3) 

and in the case of two sublattices has two values T, a (uj) 
and T, B (u>). As a result, the DMFT equations become 
matrix equations for the CDW phase. Hence, we can 
write the solution of the Dyson equation (in momentum 
space) in a matrix form 

G fc M = [z(w)-t fc ]-\ (4) 

where the irreducible part z(uj) and the hopping term t& 
are represented by the following 2x2 matrices: 



with the band structure e/- satisfying = 
— t* liniD^oo YliLi cos fci/V^D- Then we can repre- 
sent the local Green's function on sublattice a 
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in terms of the local dynamical mean field X a (uj), via 
G aa {u) = 



uj + ^ - Y, a (uj) - A a (w)' 



(7) 



Finally, we close the system of equations for Y, a (uj) and 
X a (u) by finding the local Green's function from the solu- 
tion of an impurity problem in the dynamical mean field 
X a (u). For the Falicov-Kimball model such a problem 
can be solved exactly and the result is equal to 



G QQ H 



f 
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w + /^-A Q (w) uj + fi^ - U- \ a {uj)' 



(8) 



where is the average concentration of the localized 
electrons on the sublattice a. In the CDW phase, the 
total concentration of localized electrons is fixed + 
nf = const, and the difference of the concentrations on 



each sublattice An/ = 



— rij is the order parameter 
of the CDW phase and is defined from the equilibrium 
condition on the sublattice chemical potentials: p/^ — 

M? = 0. 

Numerical solutions of these equations are given in 
Ref. HE where the evolution of the DOS in the CDW- 
ordered phase is shown. At T = 0, a real gap develops of 
magnitude U with square root singularities at the band 
edges (even on the hypercubic lattice which has infinite 
tails to the DOS in the normal state). As the temper- 
ature increases, the system develops substantial subgap 
DOS which are thermally activated within the ordered 
phase. Additional plots of the DOS can be found in 
Ref. [H]. Note that the singular behavior occurs for one of 
the "inner" band edges on each sublattice, and that the 
subgap states develop very rapidly as the temperature 
rises and completely fill in the CDW gap at the critical 
temperature T c . 



III. FORMALISM FOR INELASTIC LIGHT 
SCATTERING 

The interaction of a weak external transverse electro- 
magnetic field A with an electronic system with nearest- 
neighbor hopping is described by the Hamiltonia n 23 i 24 : 



k 



(9) 



2h 2 c 2 



E E A a {~k) la . p (k + k')A (-k'), 



kk' 



where the current operator and stress tensor for itinerant 
electrons are equal to 

Ml) = E ^T^(* + 9/ 2 H(fc - 9/2) (10) 



abk 



and 



^(9) = E ^r di{k + q/2)db{k 9/2) ' 



respectively. Here t a b(k) are the components of the 2x2 
hopping matrix in Eq. ([5]). The general formula for the 
inelastic light scattering cross section 



-Pei 



fi) 



(12) 



U 



$>(fcOfl(fc/)4ej(/|M°0(g) 

a/3 



is expresse d 23 ' 24 through the square of the scattering op- 
erator 



M a0 (q) i) = (fh«Al)\i) (13) 
+ yYaiJ>(*/)IO (l \ja(-ki)\i) 



El — Ei — U)i 



El - Si +UJf J 

which contains both nonresonant and resonant contri- 
butions. Here Q — Wi — u>f and q — ki — kf are the 
transferred energy and momentum of the photons, re- 
spectively, e l( ^) is the polarization of the initial (final) 
states of the photons and £i(f) denotes the electronic en- 
ergies for the initial i and final / electronic eigenstates. 
The quantity g(q) = (he 2 /Vug) 1 / 2 is called the "scat- 
tering strength" with uj q — c\q\, and Z is the partition 
function for the electronic system. The nonresonant part 
of the scattering operator M(q) is constructed from the 
stress tensor and the resonant one is constructed from the 
square of the current operators. After substituting the 
expression for the scattering operator into the formula for 
scattering cross section, one obtains three terms in the 
response function x(g, J7): a nonresonant term; a mixed 
term; and a pure resonant term. The result is 



R(q,n) 



27rg 2 (k l )g 2 (k f ) 
1 - exp(- ( 8n) 



(14) 



where 



X (q, SI) = X n(q, SI) + XAi(q, SI) + xn(q, SI)- (15) 
■ — I 

In Ref. 3|, we have described in detail how to extract the 
components of the cross section from the appropriate cor- 
relation functions in the normal phase: we must calculate 
corresponding multi-time correlation functions for imag- 
inary Matsubara frequencies and then analytically con- 
tinue to the real axis. Inelastic light scattering examines 
charge excitations of different symmetries by employing 
polarizers on both the incident and scattered light. The 
Ai E symmetry has the full symmetry of the lattice and is 
primarily measured by taking the initial and final polar- 
izations to be e % — e* = (1, 1, 1,1,.. .). The B lg symme- 
try involves crossed polarizers: e l — (1,1,1,1,...) and 
e* = (—1, 1, —1, 1, . . .); while the £?2g symmetry is also 
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using crossed polarizers, but with the polarizers rotated 
by 45 degrees; it requires the polarization vectors to sat- 
isfy e 4 = (V2,0,V2,0,...) and e f = (0, y/2, 0, \/2, . . .)• 
(Note in previous work we used the wrong normalization 
for the B 2g polarization vectors resulting in a resonant 
response a factor of four smaller) . For Raman scattering 
(q = 0), it is easy to show that for a system with only 
nearest-neighbor hopping and in the limit of large spatial 
dimensions, the A± g sector has contributions from non- 
resonant, mixed and resonant scattering, the B\ g sector 
has contributions from nonresonant and resonant scat- 
tering only, and the B 2g sector is purely resonan t 25 ' 26 . 
These results continue to hold in the ordered phase. 



IV. MIXED AND RESONANT 
CONTRIBUTIONS TO THE SCATTERING 
RESPONSE 

Since the nonresonant contributions to Raman scatter- 
ing in the ordered phase have already been determined^, 
we focus here on the modifications needed in the ordered 
phase to calculate the mixed and resonant responses. As 
discussed above, the mixed and resonant response func- 
tions are extracted from the corresponding multi-time 
correlation functions. For the mixed one, the appropri- 
ate response function is built on the stress tensor and two 
current operators, as follows 

X7,/,i(n,72,,7s) = (T T y(n)jWfo)j®foj) . (16) 

The symbol T T is a time ordering operator ((...) = 
Tv[e-^ H ...]/Z). Here we have introduced a compact 
notation for the contraction of the stress tensor and cur- 
rent operators [Eqs. (fTUfl and (fTTj) for q = 0, fe»(/) = 0] 
with the polarization vectors, as follows: 

7 = X] e « 7a ^ e i8' 

a/3 
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(17) 



;(/) 



respectively. The next step is to perform the Fourier 
transformation from imaginary time to imaginary Mat- 
subara frequency, and, as a result, the mixed correlation 
function is represented as a sum over Matsubara frequen- 
cies of the generalized polarizations as follows: 



(18) 



Here 



introduce the shorthand notation 
n M (iw m - ivj,iu> m + ivi - ivf,iw m ) 
for the dependence on the fermionic iu) m — iTrT(2m+ 1) 
and bosonic ivi = ilnTl Matsubara frequencies. The 



ttM 
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Ji k, iw 



FIG. 1: Feynman diagrams for the generalized polarizations 
of the mixed response function. Due to the static nature of 
the irreducible charge vertex of the Falicov-Kimball model, 
we have iuj m = iw m / . 



corresponding Feynman diagrams for the generalized 
contributions to the mixed response function are shown 
in Fig. [TJ where the first and third diagrams correspond 
to the first term in Eq. (|18[) and the other two diagrams 
correspond to the second one. 

For the resonant response function, we construct the 
four-time correlation function with four current operators 
as follows 

XiJJ,i(n,r 2) T 3 ,T 4 ) = (T T j^(n)j^(r 2 )j^(r 3 )j^(T 4 )) , 

(19) 

and in the same way as for the mixed one, the resonant 
response function is expressed as a sum of the generalized 
polarizations over Matsubara frequencies 



\, .;./.<: n',.iu>. =T^2 



(20) 



x rn fl '- 

L m, 

+ n 



R,I 

m — /,m-(-i — /,m— /' 
R,II 

m,m-\-i,m-{-i — /,m— /' 



n 



R.I 



n 



RJI I 

m,m—f,m+i—f,m+i'] ' 



The corresponding Feynman diagrams for the general- 
ized contributions to the resonant response function are 
shown in Fig. [21 where we introduce additional sublat- 
tice indices a to s. Each term in Eq. (I2U1) corresponds 
to a separate line in Fig. [2] respectively. There are also 
other contributions to the four-time correlation function 
in Eq. (flT?)) , but they do not contribute to the scat- 
tering cross section (see Ref. d for details). For the 
Bi g and B 2g symmetries, the generalized polarization 

n m,L-/,m+i-/,m-/' is a sum of tne nrst tw0 diagrams 
in the first line of Fig. [5] (the bare loop and the ver- 
tical renormalization) and the generalized polarization 

R II 

^m,m+i,m+i-f,m-f conta i ns on ly tne nrst diagram in the 
third line (the bare loop) , whereas for the A\ g symmetry, 
all diagrams in the corresponding lines (the bare loop, the 
vertical renormalization, and the horizontal renormaliza- 
tion) contribute. 

The next step is to derive analytic expressions for these 
mixed and resonant generalized polarizations. There are 
two types of Feynman diagrams (for both the mixed and 
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II 1 2 3 in the CDW phase is equal to 



loi fc, rat , 



FIG. 2: Feynman diagrams for the generalized polarizations 
of the resonant response function. 



u M,b _ ]_ y- .(*) .(/) 



(21) 



x i u fc,l ^£,2 ^fc.S + '^fc.l'J'fc^ ^,3 + 2 ( - T fc, 3 

pAB^B^AB , ^BA^tBA^BA , n BA n BB n AA 
^fc.l ^/e^ u fc,3 T ^fc,! u k,2 u k,3 T ^fcj <-*fc,2 ^fe.S 



^BS^-fAA^rSA , pBBpABpAA 



(iwi + /xf - Sf )(tw 3 + /4 - Eg 4 ) + (twi + /if 

X (i«a + Md ~ ^2 ) + (tW2 + /if - Sf )(1W 3 + Md - 

+ (i«a + Md - S^)(iw 3 + 14 - ) + (iwi + Md - 
X (iwa + /if - Ef ) + (iwi + /id - )(<W 3 + /if 



Ef) 

S3 4 ) 
Ef) 



ZiF ao (Zi) 



{2l-z\){zi-z\) {zl-zi){zi-zi) 

■^3-^00(^3) 



+ 



+ 



■^2-^00(^2) 



(zl-zi)(2i-zi) 

ZfF QO (Zi) 

(zl- z\){zl - zl) ^ (zl - zim zj) 

Z 3 F ao (Z3) 



{z\-zi){zi-zi) 



Here, we use the shorthand notation for frequencies: icJi, 
iwa, and 2 w 3 -> 1, 2, and 3, with = E a eL (/) ff^ 

and 7^ = £ Q/3 e L dLdk g e ^ and defined b y 



Z(w) = J [w + i4 — Y, A {u))][u) + (J>d — E B (o;)], (22) 



where 



Foo[^(w)] = / dep(e) 



Z(w) 



(23) 



is the Hilbert transform of the noninteracting density of 
resonant contributions as well as for the nonresonant states, which satisfies p(e) = exp(-e 2 /t* 2 )/t* y/n for the 
one±&): bare loops and renormalized loops (see Figs. Q] infinite-dimensional hypercubic lattice, 
and [2|) . First we consider the bare loops and then the 
renormalized ones. The bare term for the mixed response 



The bare loop for the resonant response n^ 2 h 3 4 (Fig. [2]) 
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in the CDW phase is equal to 



and 



A.B A,B 



n 



R.b 

1,2,3,4 



1 vrVi w i (/) i (i) i (/1 

/-^i 2-^i 2-^i •>h Jk Jk Jk 



r^na /~ibc s~ian /~iyi 
X( ~ r k,l { ~ r k,2 { - r k,3 U k,3 

X [»W3 + ^-E^][iW4 + /if -Ef] 
+ + Hd - Sf ] [»W2 + ~ ^2 ] 
x [zw 3 + /if - Sf ] [iw 4 + Md ~ Ef ] ) Xl (Z 1 ,Z 2 ,Z 3 ,Z. 



(24) 



X! + - E^][w + /if - £f ]xi(-£i> ^2, ^3, ^4 



+2 X , 1 , (Z 1 ,Z 2 ,Z 3 ,24 



Here we introduce three quantities xi{Zi, Z 2 , Z 3 , Z4), 
X l 1 {Zi 1 Z 2 ,Z z ,Z ll ), and Xi C^i> ^2, ^3, Z4), which arc 
equal to 



Xi(Zi, Z 2 , Z 3 , Z4) 

-T 



(25) 



F 00 (Z 2 )/Z 2 



{zl-zl){zl~zl){zl-ziy 



(Zf 


-zl){zl-zl){z\ 


-zl) 




F^{Z 3 )jZ 3 






— z%)(z% — Z 3 )(Z A 


-ZD 




Foo [Z/O/Zi 





Xi(Zi,Z 2 , Z 3 , Z4) 

-y 



(26) 



N V " 4)^1 el)(Z$ el){Zt e%) 
ZiF OQ (Zi) 

\zl-zl){z\-zl)(zl-zl) 



Z 2 F 0O (Z 2 



(Zf 


-zD(zi~ 


ZD{z'i 


~zD 




Z 3 -Foo 


(Zs) 




izl 


-^3 2 )(^2 2 - 


zD(zl 


~zD 




Z^Fqc 


(Z4) 





{2l-zl){zl-zl){zl-zD' 



x '{(ZuZ 2 ,z 3 ,z 4 ) 
=-y. 



(27) 



ZfFoo^Zi 



(z$-zD(z 2 3 _- z\){z\-zl) 

Z^FgciZi) 

+ {zl-zD{_zl-zD{zl-zD 
1 -^1^00(^3) 

ZlFcoiZj) 

+ (z*-z2)(zi-zl)(zi-ziy 

respectively. 

The renormalized loops in the Feynman diagrams 
in Figs. Q] and [5] describe the charge screening effects 
through the reducible charge vertex, which is defined 
through the irreducible one by a Bethe-Salpeter equa- 
tion. In the DMFT approach, the irreducible charge ver- 
tex r a is local but is different for different sublattices in 
the CDW ordered phase (see Ref. [l6l ). Nevertheless, it 
has the same functional form (when expressed as a func- 
tional of the Green's function and self-energy) as in the 
normal state^Zr— and is equal to 



pa _ 1 

1 m,m+l rp Q aa 



- "mm"- m,m+l 
■ya 



(28) 



f^aa 
"m+Z 



for the Falicov-Kimball model (an explicit formula for 
other models is unknown). This expression also fol- 
lows from the partially integrated Ward identity derived 
by Janis^S,. Because in the CDW phase the irreducible 
charge vertex is local both in the lattice and sublattice 
indices, the reducible one depends on two sublattice in- 
dices and is defined by the Bethe-Salpeter equation 

(29) 

where we introduce the bare susceptibility 



\ ^ s~iab s~iba 



m-\-l ' 



(30) 



The lattice Green functions can be derived from the 
Dyson equation in Eq. (jj) and are equal to 



k 



AA _ + ^ E^ 



(31) 



rn c fc 



s~tBB _ + H — E^ 

^fc,™ — 72 ,2 
_ € k 

/~<AB _ ^BA _ e fe 



fe,m fc,m y 2 2 * 

' ( k 
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Expressions for the renormalized loops have a similar 
form for all contributions (nonresonant, mixed, and res- 
onant) and differ only in the loops attached to the left 
and right sides of the total reducible charge vertex. The 
renormalized loop for the mixed response is then equal 
to 



n i 



M,r 
2,3 



[ xtj(iui,iu 2 ,iu)3) xfj(iui,iu)2,iu3) ] (32) 



x T 



fiAA f^AB 
1 1,3 1 1,3 



1 1,3 



BB 



x|(«wi,iu> 3 ) 



where we introduce the quantities 



xf ? -(iwi,KJ 2 ,iw 3 ) 



N 2^ J k Jk 



(33) 



\r<AA r iBA 
iAB riAA 



BA 



BA 



(~il\B riJ\J\ /~ti>J\ , 



i{uj\ + cj 2 + k>s) 



k 

AA /-iBB /-tAA 
k,uii ^- T k,u)2 k,u)3 

AB (~iAB n AA I 
B ^B ^B v-iBl 



3/^f — Sf — Sf — S3 ] 
^2^00(^2) 



{Z\-Zl){Z\ 



zD 



ZiF oc (Z 3 



{z\-zl){zi-zl) 

[iui\ + - Ef ] [iw 2 + 
F«,(Zi)/Zi 



/4 



S2 4 ] [tw 3 

F 0O (Z 2 )/Z 2 



B 



-Z\){Zl-Z\) 
F oc {Z 3 )/Z 3 



{z\-zl){zi-zi) 



{z\-zl){zl 



zl) 



and 



xfj(wi,iuj 2 ,iu}3) 



1 Vl (4) 7 (/) 



(34) 



[pBA pBA /iBB 
L < - T fc,wi LT fc,W2 LT fc,w 3 
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Ac 
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(36) 



Z 3 



to the right of the charge vertex. Now we can find the 
exact expression for the vertex corrections defined by 
Eq. (|32|) with the following form 



u l,2,3 — 



1 



Ai, 3 



(37) 



x xfj (twi, <w 2l J w 3 )Tit 3 x 1 1 f TTf 3 xf (iwi, z^ 3 ) 

+ x£(*wi ) tw 2> iw 3 ) (1 - Trf 3 xff) TT^axf (twi.iws) 

+ x^(«wi, twa.iwa) (l - TI^x^) TTf 3 xf (iwi.iwa) 

1 J3 / ■ ■ ■ \rr^ T ^B BArriT-.A A/- ■ \ 

+ Xjj{lUl,lU 2 ,lL) 3 )Tr i 3 Xl,3 TT 1,3X7 («Wi, «W 3 ) 



where 



Ai, 3 = (1 - rr^x^) (1 - rrf 3 xff 

JA 
>3 



(38) 



Tr^x^frrf 3 xf A 



For the resonant response function, the renormalized 
loops in Feynman diagrams are defined in the same way 
as the mixed one and in a compact form we have 



n 



R,r 
1,2,3,4 



1 



n,3 



(39) 



A/- ■ ■ Srrrr^A ABrriT^B Bf- ■ • \ 

X Xj3\U}t,W2, lU} 3 )TT 13 Xl,3 Tr i, 3 Xjj{iu 3 , lUl) 

+ xfj{iui,iu 2 ,iu 3 ) (1 - rrf 3 xff) TTf- >a Xjj(iuia,iWi,iwi) 

+ xf 3 (iui,iuJ2,iuJ 3 ) (l - TT? >3 Xi,i) TT i, 3 xf 3 (iu 3 , iuj 4 , iwi) 

+ xfj (wi, «W2, i^ 3 )Trf 3 xf ^rr^x/j (zw 3 , iw 4 , iwi) 

Now we have the same quantities xfj (M, iw 2 , *w 3 ) to the 
left and to the right of the charge vertex. For nonresonant 
scattering, the renormalized contributions have the same 
form with Xjj(i^>ii * w 2, «w 3 ) replaced by xti^iii^) ( see 
Ref.[li). 

The total expression for the mixed generalized polar- 
ization is finally obtained as the sum of both the bare 
and renormalized contributions: 



1,2,3 



n M,b 

1,2,3 



n M,r 
1,2,3 



(40) 



on the imaginary axis. Now we have to perform an an- 
alytic continuation to the real axis. First we replace 
the sum over Matsubara frequencies by an integral over 
the real axis. Next we analytically continue Matsub- 
ara frequencies to the real axis in the following order: 



first ? 



IVf = 



IV4 — IV 



f 



UJ 



± 



fl ± i0+ followed by 
± i0 + , and finally 
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Auj = uj[ - u)i = oj'f - Uf ->• in Eq. (|T8)1. Then the 
mixed response function is expressed directly in terms of 
the generalized polarizations as 



the Bi g and Bi % channels, the generalized polarizations 
YiRJi con tain only the bare loop contributions (the first 
diagrams in the last two lines of Fig. 
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(43) 
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On the other hand, the generalized polarization II ' con- 
tains both the bare and vertically renormalized contri- 
butions (the first two diagrams in the first two lines of 
Fig. [2]) in the B\ g and i?2 g symmetry channels: 



RJ 

Bi g , 1,2, 3, 4 
RJ 

B 2g , 1,2,3,4 



11 1.2,3,4 ' ii l,2,3,4 ) 



n 



2,3,4 
RJ 

Bi„, 1,2,3,4 • 



(44) 



In the A\ g channel, all diagrams in Fig. [2] contribute, 
hence 



n 



where f{oj) = 1 /[cxp(f3uj) + 1] is the Fermi-Dirac dis- 
tribution function. Since the imaginary- axis form of the 
response is expressed as a functional of the Green's func- 
tions and self-energies, one simply replaces the appropri- 
ate Matsubara frequency arguments by the real frequen- 
cies, according to the different terms listed above. This is 
a tedious, but straightforward exercise to yield the final 
formulas, which are too cumbersome to include here. 

For the resonant response function, an analytical con- 
tinuation onto the real axis is more complicated, but the 
general approach remains the same and final expression 
is the following: 
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(45) 



R,r 
1,2,3,4 



+ n 



R,r 

2,3,4,1- 



l 



(27ri) 2 

x \ lim [ 
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U RJ (uo-iO 
jRJ { , , ,. i0 

iO + 



duj \f(u) - f(oj + ft)} 



(42) 



It should be noted that some renormalized terms in 
Eq. (|42[) contain nominal divergences in the limit Auj — > 
[connected with vanishing determinants in Eq. (|38l) which 
are found in the denominators of Eq. but the con- 

tribution of these terms to the response is actually finite. 
In the case of the uniform phase of the Falicov-Kimball 
model, their contributions were calculated analytically 
using l'Hopital's rule^, but in the case of the CDW phase, 
the expressions are more cumbersome, so we calculate the 
limit Alo — > numerically. 



V. NUMERICAL RESULTS 
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Now that all of the formal developments are complete, 
we are ready to discuss the numerical results found by cal- 
culating the total Raman response function for different 
symmetry channels and different interaction strengths as 
functions of T within the ordered phase. We shall con- 
sider two cases: the case of a weakly scattering metal in 
the normal state (U = 0.5, T c — 0.0336) and the case 
of a strongly correlated insulator in the normal state 
(U = 2.5, T c — 0.0724); both cases are insulators at 
zero temperature due to the CDW order. In previous 
workiii, we have calculated the temperature evolution of 
the single particle DOS in the CDW phase of the Falicov- 
i0 + )] Kimball model. Here we present figures of the DOS for 
the temperatures that we calculate the Raman response 
(all temperatures are below T c ): T = 0.02 for the case of 
fl - i0 + ,oj - uj f + £0+) U = 0.5 (Fig. [3]) and T = 0.06 for the case of U = 2.5 
i0+ oj - uj f - i0 + , uj + n + M+, uj + uj, - i0+) ( Fi S-H> respectively. One common feature of the CDW- 

f rdered DOS is the presence of a sharp inverse square- 
oot-like feature at U/2 for sublattice A and -U/2 for 
sublattice B, which frame the gap as T —} 0. There 
Now we can specify the different contributions to the also are bands of subgap states with a maximum DOS 
resonant response in the different symmetry channels. In at ±E/2; we have E « 0.18 for U — 0.5 and E w 1.7 
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FIG. 3: (Color online.) Conduction electron DOS at T = 0.02 
for U = 0.5. The solid black curve is the total DOS, while the 
dashed red line is for the A sublattice and the dot-dashed blue 
line is for the B sublattice. Note how there is a divergence at 
the band edge on each sublattice which develops as T — > 0, 
and that the subgap states disappear as T — > 0. Finally, we 
have marked the locations of the band edge at ±17/2, and of 
the peak of the subgap states at ±B/2. 



3 4 
< 





-U/2 


U/2 






i 




A 












B ; 












total ; 




J 


-E/2 


E/2 










i 


1 



















-3 -2-10123 
Frequency u)/t 

FIG. 4: (Color online.) Conduction electron DOS at T = 0.06 
for U — 2.5. The solid black curve is the total DOS, while the 
dashed red line is for the A sublattice and the dot-dashed blue 
line is for the B sublattice. Note how there is a divergence at 
the band edge on each sublattice which develops as T — > 0, 
and that the subgap states disappear as T — > 0. Finally, we 
have marked the locations of the band edge at ±17/2, and of 
the peak of the subgap states at ±E/2. The DOS has upper 
and lower Mott shoulders for the strongly correlated system 
(indicated by unlabeled arrows). 



for U = 2.5. These subgap states originate from thermal 
excitations of the CDW order and they vanish at zero 
temperature. In addition, for the case of the strongly 
correlated insulator U = 2.5, one can observe in Fig. ^ad- 
ditional shoulders at ±1.82, which are resulting from the 
upper and lower Hubbard bands of the high temperature 



normal state Mott insulator. At zero temperature the 
states below zero energy are filled and the states above it 
are empty. At finite temperature, due to thermal occupa- 
tion, there are some empty states below the chemical po- 
tential and some occupied ones above. These thermally 
activated states give contributions to the optical conduc- 
tivity and to the nonresonant Raman scattering, creating 
different peaks in those function o 15 ' 16 . A large main peak 
at U corresponds to single-particle transitions from the 
lower occupied to the upper empty bands of the CDW, 
which are separated by a gap of width U. This peak will 
become more enhanced as T — > 0. Peaks also occur at 
(U +E)/2, which correspond to single-particle transitions 
from the lower occupied CDW band at —U/2 to the upper 
empty subgap states at E/2 and from the lower occupied 
subgap states at —E/2 to the upper empty CDW band 
at U/2. In addition, we see peaks at E corresponding 
to transitions from the lower occupied subgap states at 
—E/2 to the upper empty subgap states at E/2. The in- 
tensity of these peaks will decrease as T is lowered, since 
the subgap states will lose spectral weight, and eventually 
vanish. Finally, there is an additional peak at (U — E)/2, 
which corresponds to transitions between the almost fully 
occupied lower CDW band at —U/2 and the lower sub- 
gap states at —E/2 and between the almost empty upper 
subgap states at E/2 and the upper CDW band at U/2. 
The intensity of this peak will also shrink as T is low- 
ered. We anticipate all of this structure will also to be 
seen in the total electronic Raman scattering, but the 
details of the temperature dependence, or of the reso- 
nant effects are difficult to guess without performing the 
calculations. We do see, however, that we have a wide 
number of different "gap edges" where one would expect 
large resonant effects. The largest should occur when the 
photon energy is equal to U, but we should also see them 
at (U±E)/2 and E. 

Analysis of the expression in Eq. (l4"2l gives that, in 
addition to the nonresonant peaks at fi = U, (U + E)/2, 
E, and (U — E)/2, there can also exist peaks which orig- 
inate from two particle transitions, i.e. ft = (3E — U)/2, 
(U + E)/2, E, U-E, and (U-E)/2, some of which coin- 
cide with single particle transition energies. In addition, 
there can be strong resonant enhancement when either 
LOi or u>f approach these energies. 

In Fig. [5J we plot the total Raman response at T = 0.02 
for the Bi g symmetry channel with U = 0.5 for different 
energies of the incident photons. This case corresponds 
to a moderately correlated metal in the high-temperature 
phase, with a CDW gap of size 0.5. The total Raman re- 
sponse function for the Bi g symmetry contains two con- 
tributions: the nonresonant contribution (dashed line), 
which is the only contribution at very high photon en- 
ergies u>i — > oo, and the resonant contribution, which is 
also the total (resonant) response for the B^g symme- 
try. For small values of a>j, we observe only a continu- 
ous enhancement of the spectra until Wj is large enough 
to create excitations across the smallest subgaps in the 
thermally excited DOS. The uj; L = 0.3 and uii = 0.346 
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FIG. 5: (Color online.) Total Raman spectra for Bi s sym- 
metry for different values of the incident photon frequency 
and for different vertical scales in the different panels with 
T = 0.02 for U — 0.5. Colors are used for the different in- 
cident frequencies, which can also be read off by examining 
the location of the unphysical divergence when £7 — > uii on 
the hypercubic lattice. The nonresonant response (the case 
of uji = oo) is also shown with a dashed line. 



FIG. 6: (Color online.) Total Raman spectra for Bi s sym- 
metry for different values of the incident photon frequency 
and for different vertical scales in the different panels with 
T — 0.06 for U — 2.5. Colors are used for the different in- 
cident frequencies, which can also be read off by examining 
the location of the unphysical divergence when Q, — > uii on 
the hypercubic lattice. The nonresonant response (the case 
of cji — oo) is also shown with a dashed line. 



curves correspond to the initial transition of the electron 
from the lower CDW band to the upper subgap states 
with a further transition to the lower subgap states with 
an energy loss around Q, ~ (U — E)/2 = 0.16 and from 
the lower subgap states to the upper CDW band with a 
further transition to the upper subgap states with an en- 



ergy loss around fl ~ E = 0.18 [see panel (a) for details]. 
In addition, there is a peak which corresponds to the 
two particle excitations around f2 ~ (3E — U)/2 = 0.02. 
When the energy of the incident photons is tuned out 
resonance with these subgap states (e.g. uii = 0.4), the 
intensity of the peaks rapidly decreases until we approach 
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FIG. 7: (Color online.) Resonant profiles in a semilog plot for 
different values of the transferred photon frequency at T — 
0.02 for U = 0.5. Different colors denote different transferred 
frequency fl (which can also be found from the unphysical 
divergence at oJi — > Q). Note how similar the response is for 
different symmetry channels. 



the next resonance at u>i = U, which corresponds to the 
initial transitions from the lower to upper CDW bands, 
with further transitions to all states below. In this case, 
we observe the largest resonant enhancement [see panel 
(c)] of more than a factor of 1000. Note that we also have 
"joint" resonance effects, as there are multiple peaks res- 
onanting with this incident photon energy, but the res- 
onance rapidly decreases and becomes small again once 



■ i i i i 



■ i i i i i i i i i i i i i i i i i i i i i i i 




a 



12 3 4 

Incident photon energy 0)./t 

FIG. 8: (Color online.) Resonant profiles in a semilog plot for 
different values of the transferred photon frequency at T = 
0.06 for U = 2.5. Different colors denote different transferred 
frequency Q (which can also be found from the unphysical 
divergence at u>i — > Q). 



OJi reaches about 0.7 [see panel (b)]. Increasing the inci- 
dent photon energy further leads to a continuous decrease 
of the resonant response without any significant change 
in its shape; the high energy peak simply moves to the 
higher frequencies and the response settles into the non- 
resonant one. Note that every curve shows a large peak 
in the limit where f2 — > u>i. This peak is an artefact of the 
infinite-dimensional limit and the hypercubic lattice, and 
is not expected to be seen in any real material system. 

Similar behaviour is observed for the case of a strongly 
correlated insulator (in the normal state) at U = 2.5 and 
T = 0.06 in Fig. [5] The main differences with the previ- 
ous case are connected with two points. First, the gap is 
larger and the subgap states are wider separated. Hence, 
the response is very small in the low-energy part of the 
spectrum and for low initial photon frequencies. Second, 
the single particle excitation energies are quite different. 
For the case of U = 0.5, the energies of the single par- 
ticle excitations (U — E)/2 = 0.16 and E = 0.18 are 
close to each other and the corresponding peaks of the 
response functions effectively merge. Now these peaks at 
(U — E)/2 = 0.4 and E = 1.7 are well separated and 
can be distinguished in the spectrum. In addition, as 
was seen for the A\ g total Raman response in the normal 
state of the Falicov-Kimball model^, the mixed contribu- 
tion becomes large and negative for large enough values 
of the transfered frequency f2 and can completely cancel 
the resonant contribution when one is in the Mott insu- 
lator phased. Moreover, for some values of il the sum 
of the mixed and resonant contributions is negative and 
the total Raman response for the A\ g symmetry becomes 
smaller than the nonresonant one (not shown here). 

Another important feature to examine in the total Ra- 
man response is the resonant profile of the response, 
which is a cut through the spectra with a fixed value 
of the transferred energy f2 while varying the incident 
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photon frequency Wj. In Figs. [7] and[8j we plot the to- 
tal Raman response functions for different symmetries at 
various (fixed) transferred frequencies ft as a function of 
the incident photon frequency u>i. 

In the case of U = 0.5, for small values of the trans- 
ferred frequency Q — 0.1, we observe a wide peak cen- 
tered around coi ~ 0.3, which correspond to the joint 
resonance when uj; L ~ (U + E)/2 is tuned to the single 
particle transitions from the lower CDW band to the up- 
per subgap states and u)t = u)i — Q, ^ E is tuned to 
transitions between the lower and upper subgap states. 
Another sharp peak at uji = 0.6 corresponds to transi- 
tions with the scattered frequency wj — U. For larger 
values of the transferred frequency i7, the resonant pro- 
files become more complicated and dramatically change 
as the transferred frequency is increased. This compli- 
cated behavior is caused by the requirement to satisfy 
the resonance conditions when the frequencies u>i, Uf, 
and — LUi ~ ujf must be tuned to the available single 
particle transitions. Due to this constraint not all of the 
main resonances are seen, like the one at u)i = U. But 
for the large values of the transferred frequency O > U, 
when only transitions between the lower and upper CDW 
bands are involved, the shape of the resonant profiles 
changes smoothly and slowly. 

For large values of U = 2.5, when the peaks of the 
single particle DOS (as well as the energies of the single 
particle transitions) are well separated, the resonant pro- 
files display much more complicated behavior (see Fig. [51 
where we show just one symmetry channel, since all chan- 
nels are very similar on the log scale) . The overall profiles 
are significantly enhanced when the transferred frequency 
is larger than about 0.8. The profiles also illustrate peaks 
which change shape dramatically as f2 is changed. Such 
behaviour is similar to what was seen for the resonant 
profiles in the normal stated. 

VI. CONCLUSIONS 

In this work, we have shown how one can solve for the 
exact total electronic Raman response of a CDW insu- 
lator that is formed via a nesting instability. Since the 



DOS reconstructs significantly below T c , we also see a sig- 
nificant change in the Raman response as a function of 
T. The exact solution is made possible for the Falicov- 
Kimball model in the infinite-dimensional limit, where 
DMFT is exact. We use the Falicov-Kimball model be- 
cause the charge vertex is known exactly for this model. 

Our main results are that there are a large number 
of strong resonances associated with all of the different 
peaks in the ordered-phase DOS, which has significant 
subgap states at low T. The strongest resonance occurs 
between the states separated by U corresponding to the 
T = gap. Since most CDW systems have gaps less than 
an electron volt, this resonance would not normally be 
able to be seen with optical light. If the incident photon 
frequency is larger than the gap, we can, nevertheless, 
see some joint resonances, where lower-energy peaks res- 
onate, similar to what was seen in in previous normal 
state calculations. In any case, we feel these results in- 
dicate that there should be very interesting Raman scat- 
tering structures seen in experiment when one examines 
resonant effects in materials where the ordering yields a 
divergence in the single-particle DOS at T = 0, such as 
the CDW case we examined here. Hopefully, these kinds 
of experiments will be undertaken soon. 

More interesting is the case when the incident photon 
energy can be on the order of the CDW gap. To do this, 
we need to find materials with larger gaps than most 
currently known CDW systems. Perhaps these kinds of 
materials can be found in the future and the experiments 
we envision carried out on them as well. 

In any case, what is clear is that resonant effects to 
electronic Raman scattering in ordered systems can yield 
a wide range of interesting results, even if a microscopic 
description of the physical behavior is challenging. 
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